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(57) ABSTRACT 

The present invention provides methods and apparatus for 
determining a dynamical property of the systemic or pulmo- 
nary arterial tree using long time scale information, i.e., infor- 
mation obtained from measurements over time scales greater 
than a single cardiac cycle. In one aspect, the invention pro- 
vides a method and apparatus for monitoring cardiac output 
(CO) from a single blood pressure signal measurement 
obtained at any site in the systemic or pulmonary arterial tree 
or from any related measurement including, for example, 
fingertip photoplethysmography. 

According to the method the time constant of the arterial tree, 
defined to be the product of the total peripheral resistance 
(TPR) and the nearly constant arterial compliance, is deter- 
mined by analyzing the long time scale variations (greater 
than a single cardiac cycle) in any of these blood pressure 
signals. Then, according to Ohm’s law, a value proportional 
to CO may be determined from the ratio of the blood pressure 
signal to the estimated time constant. The proportional CO 
values derived from this method may be calibrated to absolute 
CO, if desired, with a single, absolute measure of CO (e.g., 
thermodilution). The present invention may be applied to 
invasive radial arterial blood pressure or pulmonary arterial 
blood pressure signals which are routinely measured in inten- 
sive care units and surgical suites or to noninvasively mea- 
sured peripheral arterial blood pressure signals or related 
noninvasively measured signals in order to facilitate the clini- 
cal monitoring of CO as well as TPR. 
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METHODS AND APPARATUS FOR 
DETERMINING CARDIAC OUTPUT 

CROSS-REFERENCE TO RELATED 

APPLICATION 5 

This application claims priority to U.S. Provisional Patent 
Application 60/446,385, filed Feb. 10, 2003, which is hereby 
incorporated by reference. 

to 

GOVERNMENT SUPPORT 

This invention was made with government support under 
Contract Number NAGS-4989, awarded by the National 
Aeronautics and Space Admini stration (NASA). The govern- 1 5 
ment has certain rights in this invention. 

BACKGROUND OF THE INVENTION 

Cardiac output (CO) is defined to be the volume of blood 20 
ejected by the heart per a unit time. Since CO also represents 
the total flow of blood supplying all the tissue beds of the 
body, it is perhaps the most indicative quantity of the state of 
the heart and circulation. CO is routinely measured in inten- 
sive care units and surgical suites in order to monitor and 25 
guide therapy for critically ill patients. These patients include, 
for example, those in shock (e.g., cardiogenic, hemorrhagic, 
or septic) or heart failure and those during and after surgery 
(e.g., coronary artery bypass grafting or heart valve replace- 
ment.) 30 

An ideal CO measurement technique would be simple to 
perform, inexpensive, noninvasive or minimally invasive and 
safe, and very accurate. However, none of the conventional 
measurement techniques known in the art possess all of these 
characteristics [10]. For example, the thermodilution tech- 35 
nique, which is currently employed in most intensive care 
units and surgical suites, involves injecting cold saline into 
the right atrium and measuring the temperature downstream 
in the pulmonary artery. The average CO over the measure- 
ment period may than be computed from conservation of 40 
mass laws. Although the technique is relatively simple and 
inexpensive, it requires an invasive right heart catheterization 
whose safety is questionable [8,38] and is not very accurate 
due to the many assumptions upon which it is based (e.g., no 
saline recirculation and thorough blood mixing) [ 1 0,27] . The 45 
most accurate, conventional technique for measuring CO 
involves surgically implanting a flow probe, either electro- 
magnetic or ultrasonic, directly on the aorta. Although this 
technique also provides a continuous measurement of CO, it 
requires a high risk thoracotomy which is rarely performed in 50 
practice. Moreover, the accuracy of the aortic flow probe is 
highly dependent on vessel preparation and may only be 
accurate to within about 15-20 percent [10,27]. 

Although the development of an ideal CO measurement 
technique has proven to be difficult, several ideal, or near 55 
ideal, techniques are currently available for the continuous 
measurement of peripheral arterial blood pressure such as 
Finapres technology [23] and arterial tonometry [25]. Previ- 
ous investigators have therefore sought techniques to monitor 
CO from peripheral arterial blood pressure signals. The most 60 
popular techniques in the art are the so-called pulse contour 
methods that assume the arterial tree to be well represented by 
a parallel combination of a capacitor and resistor thereby 
accounting for the compliance of the large arteries (AC) and 
the total peripheral resistance (TPR) of the small arteries. If 65 
the instantaneous CO supplied by the heart is represented as 
a current source, then the simple model of the heart and 
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arterial tree in FIG. 1A results. Most of these types of pulse 
contour methods are specifically based on mathematical for- 
mulas which are derived by making simplifying assumptions 
and approximations to this model (see, for example, [7, 
19-21, 26, 40, 43, 44]). These methods have generally failed 
to yield good correlation between CO determined from analy- 
sis of ABP signals and directly measured CO over a wide 
range of physiologic conditions [39, 42]. 

Bourgeois et al. [3,4] did successfully demonstrate that 
their pulse contour method when applied to an ABP signal 
measured centrally in the aorta could yield a quantity which 
varied linearly with directly measured CO (electromagnetic 
aortic flow probe) over a wide range of physiologic condi- 
tions. Their method, which makes no simplifying assump- 
tions or approximations to the model of FIG. 1A, may be 
explained as follows. During the diastolic period of each 
cardiac cycle, the heart is filling and not supplying blood to 
the arterial tree (see FIG. IB). Thus, according to the model, 
ABP decays with a time constant x D equal to the product of 
TPR and AC during each diastolic period (see FIG. IB). Since 
AC is essentially constant over a wide pressure range and on 
the time scale of days [4, 1 8,37], CO could then be monitored 
to within a constant scale factor (equal to the reciprocal of 
AC) by dividing the ABP by x D . 

Bourgeois et al. specifically demonstrated that beat -to-beat 
CO may be monitored from x D and the governing differential 
equation of the model of FIG. 1 A. Thus, a key step of the pulse 
contour method of Bourgeois et al. is to fit an exponential to 
the diastolic decay portion of an ABP wavelet in order to 
measure x D . Osborn et al. [34] introduced essentially the 
same method prior to Bourgeois et al. but their experimental 
validation was not as complete or compelling. 

Bourgeois et al. were able to validate their pulse contour 
method with respect to a canine ABP signal measured cen- 
trally in the aorta, because the diastolic portion of such a 
signal usually resembles an exponential decay (see FIG. 2A). 
These investigators specifically identified the position in the 
aorta at the level of the dorsal insertion of the diaphragm as 
the optimal site for observing an exponential diastolic decay. 
However, Bourgeois et al. acknowledged that central ABP is 
rarely obtained clinically because of the difficulty in inserting 
a catheter retrogradely via a peripheral arterial blood vessel 
and the risk of blood clot formation and embolization. More- 
over, they recognized that, in peripheral ABP signals which 
are routinely made available in intensive care units and sur- 
gical suites usually via a more simple and safe radial artery 
catheterization, an exponential diastolic decay is usually not 
apparent (see FIG. 2B). The method of Bourgeois et al. there- 
fore cannot generally be applied to readily available periph- 
eral ABP signals. In fact, its application to central ABP sig- 
nals may be somewhat limited, as Cundick et al. [9] reported 
that they could not identify an optimal location in the human 
aorta in which the diastolic portion of the ABP signal 
appeared as a pure exponential decay. 

Other pulse contour methods that are based on more com- 
plex representations of the arterial tree have also been devel- 
oped (see, for example, [5, 1 1, 12, 16, 24, 31-33, 50]). How- 
ever, these techniques required the analysis of one, or even 
two, central ABP signals. Thus, their clinical utility is also 
severely limited. 

Several techniques have more recently been introduced in 
an attempt to monitor CO from ABP signals measured periph- 
erally. Techniques based on an adaptive aorta model which 
require ABP signals measured at two peripheral sites — the 
carotid artery and the femoral artery — have been developed 
[36, 46]. However, catheters are usually not placed for pro- 
longed periods of time at either of these sites in intensive care 
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units or surgical suites due to issues of safety. Another previ- 
ous technique is based on an empirically-derived formula 
which involves the calculation of the derivative of the ABP 
signal [14]. However, in order to mitigate the corruptive 
effects of wave reflections on the derivative calculation, this 5 
technique also requires two peripheral ABP measurements, 
one of which is obtained from the femoral artery. Other tech- 
niques based on a learning approach have been previously 
proposed [6, 15, 30]. However, these techniques require 
extremely large sets of training data consisting of simulta- 10 
neous measurements of CO and ABP signals obtained over 
the entire range of physiologic conditions. Moreover, the 
success of these techniques was only demonstrated with cen- 
tral ABP signals or only over a narrow physiologic range. 
Finally, Wesseling et al. [ 1, 48, 49] and Linton and Linton [28] 15 
have recently proposed model-based techniques which 
require only the analysis of a single radial artery pressure 
signal. However, Linton and Linton showed that their tech- 
nique was reasonably accurate only over a narrow range of 
physiologic conditions, and several previous studies have 20 
demonstrated the inadequacy of the method of Wesseling et 
al. (see, for example, [13, 22]). 

It is evident that there remains a need in the art for methods 
and apparatus for determining CO reliably and accurately 
using information obtained from the arterial blood pressure 25 
signal. In particular, there remains a need in the art for meth- 
ods and apparatus for determining CO reliably and accurately 
using information obtained from the peripheral arterial blood 
pressure signal. 

30 

SUMMARY OF THE INVENTION 

The present invention addresses this need, among others. 

In one aspect, the invention provides for the measuring of a 
physiologic signal indicative of cardiovascular system activ- 35 
ity, e.g., an arterial blood pressure (ABP) signal over a plu- 
rality of cardiac cycles. From analysis of the arterial blood 
pressure signal or other physiological signal the times of the 
cardiac contractions are identified. Then the relationship 
between the times of the cardiac contractions and the physi- 40 
ologic signal over a plurality of cardiac cycles is mathemati- 
cally analyzed. From this analysis a mathematical relation- 
ship between the occurrence of the cardiac contraction at a 
certain time and the subsequent time evolution of the arterial 
blood pressure over a time period greater than one cardiac 45 
cycle is obtained. This mathematical relationship is then used 
to determine a dynamical property of the system. For 
example, in one preferred embodiment the mathematical rela- 
tionship is the impulse response function between the occur- 
rence of cardiac contractions and the arterial blood pressure 50 
signal. In one preferred embodiment the dynamical property 
is the time constant which describes the decay of the impulse 
response function over a time interval, e.g., a time interval 
such as between 2 and 4 seconds following the maximum 
height of the exponential. In another preferred embodiment 55 
the dynamical property is the impulse response function 
itself. 

In particular, in a preferred embodiment the invention pro- 
vides a method and apparatus for monitoring CO by analyz- 
ing the long time scale variations (greater than a cardiac 60 
cycle) in a single ABP signal, which may be obtained at any 
site in the systemic or pulmonary arterial tree. The present 
invention determines TD through the analysis of long time 
intervals (60-90 second intervals in a preferred embodiment) 
of an ABP signal according to the following three steps (see 65 
FIG. 3). The first step involves constructing an impulse train 
signal, x(t), in which each impulse is located at the start of a 
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cardiac contraction and has an area equal to the ensuing 
arterial pulse pressure (or an arbitrary constant value in 
another preferred embodiment). The constructed signal 
therefore approximately reflects cardiac contractions in terms 
of timing and output (or just timing, in those embodiments in 
which the area of the impulses is an arbitrary constant). 

The second step deals with determining the dynamical 
properties of the arterial tree through the characterization of 
the relationship between the cardiac contractions and the 
ABP signal. This is achieved by estimating an impulse 
response function (h(t)) which when convolved with x(t) best 
fits the ABP signal (y(t)), according to any of a number of 
criteria. The estimated h(t) represents a normalized ABP 
response to a single cardiac contraction. The final step 
involves fitting an exponential to the tail end of the diastolic 
decay portion of the estimated h(t) in which the faster wave 
reflections have vanished in order to determine x D . Accurate 
determination of x D is achieved by virtue of h(t) coupling the 
long time scale variations in x(t) to y(t). 

The present invention, which in a preferred embodiment 
includes an analog-to-digital converter, a buffer system, a 
signal processing unit, and a display and alarm system (see 
FIG. 4 ), may thus be utilized to monitor CO and/or TPR, to 
within constant scale factors from any measured ABP signal 
despite the presence of wave reflections. Note that absolute 
CO and/or TPR of the systemic or pulmonary arterial tree 
(depending on the signal measurement site) may also be 
determined, if desired, by calibration with a single, absolute 
measure of CO such as a thermodilution measurement. 

The present invention may be employed in intensive care 
units and surgical suites in which invasive radial ABP signals 
are routinely available and pulmonary ABP signals are some- 
times available. The present invention may also be applied to 
noninvasively measured peripheral ABP signals (e.g., 
Finapres technology, arterial tonometry) or noninvasively 
measured signals related to ABP (e.g., fingertip photoplethys- 
mography [45], ear densitography [17]). Thus, the present 
invention could easily be employed in primary care settings, 
emergency rooms, and regular hospital beds in order to facili- 
tate the evaluation of the patient’s heart and circulatory state. 

This application refers to various patents and publications. 
The contents of all of these are incorporated by reference. 

BRIEF DESCRIPTION OF THE DRAWING 

FIG. 1 A is a simple, lumped model of the heart and arterial 
tree (either systemic or pulmonary). CO is cardiac output; 
ABP, arterial blood pressure; TPR, total peripheral resistance; 
and AC, arterial compliance. 

FIG. IB is a plot of a realistic, continuous CO signal and a 
plot of the ABP signal response to the CO signal determined 
according to the simple model of FIG. 1A. The time constant 
x D governs the dynamical behavior of the simple model. 

FIG. 2A is an example of a central ABP signal measured in 
an experimental sheep preparation. 

FIG. 2B is an example of a peripheral ABP signal measured 
in the brachial artery of an experimental dog preparation. 

FIG. 3 is an illustration of how one embodiment of the 
present invention determines x D from a peripheral ABP sig- 
nal. PP is arterial pulse pressure; R, the time of the R-wave in 
a surface electrocardiogram; x(t), an impulse train signal 
representing cardiac contractions; y(t), an ABP signal; and 
h(t), the impulse response which when convolved with x(t) 
best fits y(t). 

FIG. 4 is a block diagram of one embodiment of the present 
invention. 
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FIG. 5 is the result of a pilot experimental rabbit study 
which demonstrates the correspondence between relative CO 
changes determined from a femoral ABP signal according to 
the present invention and an aortic flow probe. CO was altered 
by various interventions (e.g., pacing, inferior vena cava bal- 5 
loon occlusion, and nitroglycerine). The solid line is the line 
of identity, and the dashed line is the line that best fits the data 
points. 

FIG. 6 shows example excerpts of measured (a) central, (b) 
“radial”, and (c) femoral arterial blood pressure (ABP) sig- 10 
nals from an experiment performed in swine. 

FIG. 7 presents results from animal 1 in an experiment in 
which the method of the invention was tested in swine. The 
top panel illustrates the aortic flow probe cardiac output (CO) 
trend (solid), estimated and calibrated CO trends from the 15 
“radial” (dash) and femoral (dash-dot) arterial blood pressure 
(ABP) signals, and the intervention and duration (underline). 
The second panel depicts x D trends estimated from the 
“radial” (dash) and femoral (dash-dot) ABP signals. The bot- 
tom two panels illustrate the mean ABP (MAP) and heart rate 20 
(HR) trends. 

FIG. 8 present results from animal 5 in an experiment in 
which the method and apparatus of the invention were tested 
in swine. Panels correspond to those in FIG. 7. 

25 

DETAILED DESCRIPTION OF CERTAIN 
PREFERRED EMBODIMENTS OF THE 
INVENTION 

The present invention encompasses the recognition that 30 
there is significant information present in the ABP signal 
when measured over long time scales (greater than a single 
cardiac cycle), and that this information, referred to herein as 
“long time scale information”, may be used to facilitate deter- 
mination of dynamical properties of the systemic or pulmo- 35 
nary arterial tree. 

Long time scale information incorporates information 
reflecting the variability of the ABP signal between beats in 
addition to, or instead of, reflecting only information about 
the ABP signal within single beats. Although long time scale 40 
information is acquired by measuring the ABP signal over a 
period greater than a single cardiac cycle (typically seconds to 
minutes), not all quantities derived from measurements per- 
formed over time periods greater than a single cardiac cycle 
include long time scale information. For example, it is com- 45 
mon to measure ABP over a plurality of cardiac cycles and 
average the measured amplitudes to obtain mean (average) 
ABP. However, mean ABP does not incorporate long time 
scale information, because it contains no information reflect- 
ing the variability between beats. Once an average value is 50 
found (e.g., by integrating the ABP signal over a time interval 
and dividing by the length of the time interval), information 
reflecting differences in length and/or amplitude of the ABP 
signal is lost. Thus the mean ABP does not incorporate long 
time scale information. 55 

Capturing the information present over long time scales 
offers a number of advantages for the determination of 
dynamical properties of the systemic or pulmonary arterial 
tree. In particular, utilizing such information provides a 
means of accurately determining such properties through 60 
measurement of the ABP signal, despite the fact that the ABP 
signal (particularly the peripheral ABP signal) is corrupted by 
wave reflections that occur at sites of impedance mismatch 
(e.g., vessel bifurcations). 

The inventors have recognized that at sufficiently long time 65 
scales, the wavelengths of the propagating waves may be 
longer than the dimension of the arterial tree. Thus the wave 
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reflections corrupt peripheral ABP signals, but only on rela- 
tively short time scales which, however, may exceed the dura- 
tion of a single cardiac cycle, leaving the signal on longer time 
scales relatively undisturbed. This implies that the lumped 
model in FIG. 1A is a valid representation of the long time 
scale dynamics of the arterial tree (either systemic or pulmo- 
nary) despite its limitations as an accurate representation of 
the arterial tree over single cardiac cycles. Moreover, the 
arterial tree is continuously being excited at time scales 
greater than a cardiac cycle by ongoing perturbations (e.g., 
breathing) and the dynamic, compensatory response of the 
regulatory system (e.g., arterial baroreflex) [2] . Thus signifi- 
cant long time scale information is generally present in ABP 
signals, and the present invention utilizes this information to 
provide methods for determining dynamical properties of the 
systemic or pulmonary arterial tree. 

The invention provides a method for determining a 
dynamical property of the systemic or arterial tree comprising 
steps of: (a) measuring a physiologic signal over a plurality of 
cardiac cycles; (b) obtaining a relationship between the tim- 
ing of a cardiac contraction and the evolution of the physi- 
ologic signal over a time period greater than that of a single 
cardiac contraction cycle by analyzing the physiologic signal 
over a plurality of cardiac cycles; and (c) using the relation- 
ship to determine the dynamical property. In general, the 
physiologic signal is a signal indicative of cardiovascular 
system activity. For example, in certain embodiments of the 
invention the physiologic signal is an arterial blood pressure 
(ABP) signal. In certain embodiments of the invention the 
physiologic signal is a signal related to the ABP signal. Such 
related signals include, but are not limited to, pressure signals 
obtained using fingertip photoplethysmography [45], ear 
densitography [17], etc. In certain embodiments of the inven- 
tion the physiologic signal is the arterial-systemic filling pres- 
sure difference (ASFPD). For purposes of description, the 
invention will be described in terms of the ABP signal, but it 
is to be understood that other physiologic signals may be 
similarly used. 

In general, a dynamical property of a system is a charac- 
teristic of the system that relates to how the system responds 
over time to a change in one or more of the parameters of the 
system. For example, a mathematical expression that relates 
future values of one or more signals generated by a system to 
past and present values of signals either sensed or generated 
by the system would constitute a dynamical property of the 
system. A characteristic time is also a dynamical property. In 
general, a characteristic time indicates the time scale of the 
temporal evolution of a function. (For example, for an expo- 
nential function given by y(t)=e“* /T , the time constant t is a 
characteristic time of the function.) 

In particular an impulse response function which enables 
one to compute the expected future values of a signal gener- 
ated by the system from past values of that signal or past 
values of other signals either generated or sensed by the 
system would constitute a dynamical property of the system. 
The characteristic time of the decay of that impulse response 
function would also be a dynamical property of the system. 
Note that the dynamical property need not fully predict how 
the system responds over time to a change in one or more of 
the parameters of the system, rather it need only be descrip- 
tive of the response. Dynamical properties of systems can 
often be estimated from analysis of signals associated with 
the system. 

The relationship between the timing of a cardiac contrac- 
tion and the evolution of the ABP signal over time may be 
obtained in a number of different ways. Determination of one 
or more dynamical properties of the systemic or pulmonary 
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arterial tree in turn allows one to obtain values for a variety of 
important parameters that characterize the cardiovascular 
system, including, but not limited to, cardiac output, total 
peripheral resistance, cardiac index, stroke volume, charac- 
teristic time constant, etc. 5 

One application of the methods described above is in mea- 
suring cardiac output (CO). The invention provides a method 
of determining cardiac output to within a constant scale factor 
comprising steps of: (a) measuring a physiologic signal over 
a plurality of cardiac contraction cycles; (b) estimating a to 
function that represents the response of the physiologic signal 
to a cardiac contraction over a time period greater than that of 
a single cardiac cycle; (c) determining a characteristic time of 
the function; (d) determining cardiac output to within a con- 
stant scale factor by dividing the magnitude of the physi- 15 
ologic signal by the characteristic time obtained in step (c). 
The invention further provides a method of determining total 
peripheral resistance to within a constant scale factor com- 
prising steps of: (a) measuring a physiologic signal over a 
plurality of cardiac contraction cycles; (b) estimating a func- 20 
tion that represents the response of the physiologic signal to a 
cardiac contraction over a time period greater than that of a 
single cardiac cycle; and (c) determining a characteristic time 
of the function, wherein total peripheral resistance is given to 
within a constant factor by the characteristic time. 25 

In certain embodiments of the invention rather than deter- 
mining the characteristic time of the function estimated in 
part (b) of the above methods, a second function that repre- 
sents the response of a different physiologic signal to a car- 
diac contraction over a time period greater than that of a single 30 
cardiac contraction is estimated, and the characteristic time of 
this second function is determined and used in step (d). 
Accordingly, the invention provides a method of determining 
cardiac output to within a constant scale factor comprising 
steps of: (a) measuring a first physiologic signal over a plu- 35 
rality of cardiac contraction cycles; (b) measuring a second 
physiologic signal over a plurality of cardiac contraction 
cycles; (c) estimating a function that represents the response 
of the second physiologic signal to a cardiac contraction over 
a time period greater than that of a single cardiac cycle; (d) 40 
determining a characteristic time of the function; and (e) 
determining cardiac output to within a constant scale factor 
by dividing the magnitude of the first physiologic signal by 
the characteristic time obtained in step (d). The methods will 
now be described in more detail. 45 

A feature common to the techniques discussed above for 
monitoring CO from continuous ABP is that the signal analy- 
sis is considered only within individual cardiac cycles. 
Because of the presence of wave reflections at these time 
scales, these techniques are limited in that they 1 ) can only be 50 
applied to central ABP signal in which the cumulative effects 
of the pulse reflections are largely attenuated; 2) necessitate 
two peripheral ABP signal measurements which are not usu- 
ally obtained clinically; 3) require a large set of training data 
obtained over a wide range of physiologic conditions, or 4) 55 
are reasonably accurate only over a limited physiologic 
range. 

In particular, methods that attempt to determine CO by 
measuring the time constant (t^) of the arterial tree, which 
could then be used to compute CO to within a constant scale 60 
factor (equal to the reciprocal of AC) have had only limited 
success for the following reasons. In peripheral ABP signals 
such as those typically available in intensive care units and 
surgical suites via a technique such as radial artery catheter- 
ization (and also peripheral ABP signals obtained via non- 65 
invasive techniques such as fingertip photoplethysmography 
or ear densitography), an exponential decay is not usually 
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apparent (see FIG. 2B). This is because, as mentioned above, 
the arterial tree is not a lumped system but actually a distrib- 
uted system with impedance mismatches throughout the sys- 
tem due to vessel tapering and bifurcations as well as changes 
in vessel caliber. Thus the diastolic portion (as well as the 
systolic portion) of peripheral ABP signals is corrupted by 
wave reflections that occur at each site of impedance mis- 
match. (Note that the complexity of these sites as well as their 
varying distances from the aorta result in reflected waves with 
large phasic differences which generally tend to mitigate the 
cumulative effects of these waves on the central ABP signal.) 
Even in the case of central APB signals it may not be possible 
to identify an optimal location in the human aorta at which the 
diastolic portion of the ABP signal appears as a pure expo- 
nential decay [9]. Similarly, efforts to locate a suitable posi- 
tion in the pulmonary arterial tree in which the diastolic 
portions of the pulmonary ABP signal appear as an exponen- 
tial decay have been unsuccessful [41]. 

Thus the presence of wave reflections has impeded devel- 
opment of accurate methods to determine the time constant of 
the arterial tree, which could then be employed to compute 
CO. As described below, the present invention overcomes this 
difficulty through the use of long time scale information to 
accurately estimate x D , from which CO may then be deter- 
mined using the following formula: 


CQ ^ {ABP){AC) 
td 

The scale factor, AC, for a particular individual may be 
determined by obtaining a single absolute measurement of 
CO, e.g., by thermodilution, and then solving for AC in the 
above formula using the estimated value for TD. Alternately, 
AC may be estimated using tables or nomograms, which are 
well known in the art and may be based on parameters such as 
age, weight, height, or particular disease status (see, e.g., U.S. 
Pat. No. 6,485,431). Additional parameters such as total 
peripheral resistance (TPR), cardiac index, stroke volume 
(SV), etc., may also be determined using well known rela- 
tionships. For example, TPR=(t Z 5 )/AC, TPR=ABP/CO, and 
SV=(CO)/heart rate. 

The CO signal estimated from the formula given above will 
reflect true CO over time periods greater than or equal to a 
single cardiac cycle, but will generately not accurately reflect 
the cardiac flow signal within a single cardiac cycle. For this 
reason, it may be desirable to average the ABP signal or the 
estimated CO signal over each cardiac cycle. Alternatively, in 
certain embodiments of the invention the ABP signal or esti- 
mated CO signal is filtered using a low-pass filter with a 
characteristic response time greater than the duration of a 
typical cardiac cycle. Another alternative is to simply average 
the ABP signal or estimated CO signal over a time period 
longer than the duration of a typical cardiac cycle (multibeat 
averaging). In general, a single-beat average approach or by 
use of a low-pass filter with a fairly short characteristic time 
would be expected to retain the most information regarding 
time variation of the CO signal as compared with simple 
averaging of the CO signal over long time periods. The time 
constant x D changes slowly in time (because over a time scale 
of tens of seconds peripheral resistance is slowly varying and 
the arterial compliance AC may be regarded as essentially 
constant), however the cardiac output itself can vary much 
more rapidly — on a beat-to-beat basis. Thus it is generally 
preferable to estimate the CO from the formula given above 
together with use of a single beat average or low -pass filtering 
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approaches — even if the constant x D is estimated from long 
epochs of data comprising perhaps tens of seconds — as 
opposed to simple averaging over time scales long compared 
to the duration of a single cardiac cycle. 

Although absolute CO may be determined using the for- 5 
mula given above, one important aspect of the invention is the 
recognition that in many circumstances it is not necessary to 
obtain a value for absolute CO. For example, in the context of 
continuous monitoring in the acute setting (e.g., in intensive 
care units), it is changes in CO rather than absolute CO that is to 
most clinically relevant. Thus determination of the propor- 
tionality constant is unnecessary, and this potential source of 
error may thus be avoided. In other words, the present inven- 
tion may be used to monitor CO (e.g., identify and quantify 
changes in CO) instead of (or in addition to) determining 15 
absolute CO. 

According to the method of the invention an analog ABP 
signal is measured invasively or noninvasively at any site in 
the systemic or pulmonary arterial tree. The analog signal is 
quantized and sampled. For example, in a preferred embodi- 20 
ment of the method the signal is quantized at 12 bits and 
sampled at 90 Hz. It is noted that these values are exemplary 
only, and one of ordinary skill in the art will readily be able to 
select other appropriate values. A signal representing cardiac 
contractions is constructed through the formation of an 25 
impulse train in which each impulse is located at the start of 
a cardiac contraction and has an area equal to the ensuing 
arterial pulse pressure, i.e., the pulse pressure that results 
from that cardiac contraction (see FIG. 3). The start of each 
cardiac contraction i s determined by detecting the onset of the 30 
upstroke of each ABP wavelet, which may be done by any of 
a variety of methods known in the art, e.g., those described in 
[51], [52], or [53]. The arterial pulse pressure for each wavelet 
is given by the difference in the maximum value of ABP and 
the value of ABP at the onset of the upstroke. 35 

Alternatively, the area of each impulse may be set to an 
arbitrary constant value. In another preferred embodiment, a 
surface electrocardiogram (ECG) is measured simulta- 
neously with the ABP signal. The two signals are quantized at 
1 2 bits and sampled at 3 60 Hz. It is noted that these values are 40 
exemplary only, and one of ordinary skill in the art will 
readily be able to select other appropriate values. The start of 
each cardiac contraction may then be established by detecting 
each R-wave of the ECG. In certain embodiments of the 
invention the constructed impulse train and ABP signal are 45 
then decimated, e.g., for purposes of noise reduction. For 
example, in one embodiment of the invention the constructed 
impulse train and ABP signal are decimated from 360 Hz to 
90 Hz. 

The single contraction ABP response (normalized by 50 
approximately the average arterial pulse pressure when each 
impulse is scaled to the ensuing arterial pulse pressure), 
which quantitatively characterizes the dynamical properties 
of the arterial tree, is determined by estimating the discrete- 
time impulse response function (h(t)) which when convolved 55 
with the 90 Hz impulse train (x(t)) “best” fits the 90 Hz ABP 
signal (y(t)) in the least squares sense (see FIG. 3). 

The impulse response function is assumed to be well rep- 
resented by an autoregressive moving average (ARMA) 6Q 
model which is given below: 


y(t) = Yj a iy(t -0 + 2 b[X ^ - ® + e 
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where e(t) is the residual error term, m and n limit the number 
of terms in the model (model order), and the set of parameters 
{a z , b z } completely defines h(t) [29]. Because the ARMA 
model is parametric, causality may be imposed, which is 
necessary for reliably estimating h(t) as x(t) and y(t) are 
related in a closed-loop fashion (i.e., y(t) also influences x(t) 
through the autonomically mediated heart rate baroreflex) 
[47]. 

For a fixed model order, the set of parameters is estimated 
from 60-90 second intervals of x(t) and y(t) through the 
least-squares minimization of the residual error term, which 
has an analytic solution [29]. The model order is determined 
by an ARMA parameter reduction algorithm that penalizes 
for the degree of model complexity [35]. Prior to estimation 
of h(t), x(t) and y(t) may be lowpass filtered in order to 
amplify the contribution of long time scale energy such that 
the least squares fit between x(t) and y(t), at these time scales, 
is prioritized. Note that any other parametric model (e.g., 
autoregressive moving average with exogenous input (AR- 
MAX) model [29]) may be employed in various embodi- 
ments of the invention to represent the structure of h(t), and 
any other minimization criterion (e.g., absolute error) may be 
utilized to find the “best” h(t). In certain embodiments of the 
invention the number of parameters in the model is selected at 
least in part based on the heart rate. 

The x D quantity is determined by finding the “best” expo- 
nential that fits h(t) over a selected time interval following the 
time of the maximum value of h(t), preferably a time interval 
in which the faster wave reflections have become minimal. 
For example, as shown in the plot of h(t) presented in FIG. 3, 
the contribution of the faster wave reflections becomes mini- 
mal at approximately 1.5-2 seconds, as evidenced by the 
decrease in fluctuations of h(t). Following this time the 
impulse response may be accurately approximated as an 
exponential. Thus in certain preferred embodiments of the 
invention the selected time interval begins approximately 1 .5 
seconds following the time of maximum h(t), more preferably 
approximately 2 seconds following the time of maximum 
h(t). For example, it has been found that a time interval of 2 to 
4 seconds following the time of maximum value of h(t) is 
suitable. Longer time intervals may also be used. Typically 
the appropriate time interval is predetermined, but in certain 
embodiments of the invention it may be selected as the mea- 
surements are being made. 

The following exponential equation is the basis of the 
resulting least squares problem where A and x D are param- 
eters to be estimated through the least squares minimization 
of w(t). 

h{t)=Ae-t/x D +w(t), 

By first log transforming h(t) over the interval of interest 
(which is always greater than zero), the optimal estimate of 
the parameters in the least squares sense may be estimated 
through an analytic linear least squares solution [4]. CO may 
then be computed to within a constant scale factor equal to 
1/AC through the ratio of the ABP signal to x D as discussed 
above. Note that TPR of the systemic or pulmonary arterial 
tree (depending on the signal measurement site) is trivially 
given, to within a constant scale factor equal to AC, by x D . 

In the embodiment of the invention described above pro- 
portional CO is determined by dividing ABP by x^. In another 
embodiment of the invention, rather than employing ABP, the 
method uses the arterial- systemic filling pressure difference 
(ASFPD). The ASFPD is determined by subtracting systemic 
venous pressure from arterial blood pressure [54]. The sys- 
temic filling pressure can either be measured, or more com- 
monly estimated. All the analyses described above can then 
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be performed on the ASFPD rather than on the ABP signal, 
including estimation of the cardiac output. In this embodi- 
ment the ASFPD is divided by the time constant to obtain a 
signal proportional to cardiac output. 

This embodiment of the invention offers a number of 
potential advantages. In the absence of cardiac contractions 
the arterial blood pressure would decay over time to the 
systemic filling pressure, thus it is not fully accurate to 
describe the decay of the ABP impulse reponse function at 
long times as an exponential that decays to zero. However, 
since in the absence of cardiac contractions the ASFPD does 
decay to zero, a description in terms of an exponential decay 
at long times is more appropriate for the impulse response 
function of the ASFPD. Furthermore, since the time constant 
of the impulse response function is estimated at long times 
when the impulse response function has already decayed 
substantially, even if the systemic filling pressure is small, 
there may be a significant difference in the decay constant 
estimated from the impulse response functions of the ABP 
compared to the ASFPD — with the analysis of the ASFPD 
providing a more accurate result. 

FIG. 4 depicts a block diagram illustrating a preferred 
embodiment of the present invention. An analog ABP signal 
is fed into an analog-to -digital converter as it is being mea- 
sured. The ABP signal may be acquired using standard meth- 
ods, such as those mentioned above. In certain embodiments 
of the invention a surface electrocardiogram is also obtained 
(ECG), e.g., via standard ECG leads. The digitized ABP 
signal is stored in a buffer system. The most recent 60-90 
second intervals of the digitized signal are transferred from 
the buffer system to a processing unit which analyzes the 
signal according to FIG. 3 in order to estimate x D . The buffer 
and processing unit may be implemented using, for example, 
any standard microcomputer running appropriate software to 
implement the mathematical operations described above. The 
software components of the invention may be coded in any 
suitable programming language and may be embodied in any 
of a ragne of computer-readable media including, but not 
limited to, floppy disks, hard disks, CDs, zip disks, DVD 
disks, etc. Outputs such as proportional CO, CO, and TPR 
may be displayed on a visual display such as a computer 
screen and/or may be printed or transmitted to a remote loca- 
tion. The ECG, and analysis thereof, may also be displayed. 
In a preferred embodiment of the system the process is con- 
tinuously repeated thereby providing the on-line monitoring 
of CO and TPR (with a delay of 30-45 seconds). Finally, in 
certain embodiments of the invention an alarm is triggered 
upon excessive decreases in CO. 

Alternatively, the AC proportionality constant may be 
computed with a single absolute measure of CO (e.g., ther- 
modilution) through the product of x D and the measured CO 
divided by ABP, as discussed above. The proportionality con- 
stant AC may then be utilized to obtain absolute measures of 
CO and TPR. Note that physiologic changes in AC due to 
disease or aging may also be monitored with multiple, simul- 
taneous measurements of absolute CO and an ABP signal. 

The present invention was evaluated in a pilot study in 
which femoral ABP and aortic flow probe CO were simulta- 
neously measured in an experimental rabbit preparation dur- 
ing various interventions known to alter CO (e.g., pacing, 
inferior vena cava balloon occlusion, and nitroglycerine). 
FIG. 5 illustrates that the resulting relative changes in CO 
with respect to baseline predicted by the present invention 
correspond to those changes determined from the aortic flow 
probe. In FIG. 5 the x-axis represents the percentage change 
in CO as measured using the aortic flow probe, while the 
y-axis represents the percentage change in CO as determined 
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using the present invention. The solid line (line of identity) 
represents results that would be achieved by a technique that 
precisely replicated the aortic flow probe measurements. The 
dashed line represents the values measured using the present 
5 invention and demonstrates a very close fit.) 

The invention was further evaluated in six experiments in 
swine, in which peripheral ABP signals and independent CO 
via an aortic flow probe were simultaneously measured over 
a wide physiologic range. Six Yorkshire swine (30-34 kg) 
to were studied under a protocol approved by the MIT Commit- 
tee on Animal Care. The animals were given intramuscular 
tiletamine-zolazepam, xylazine, and atropine prior to endot- 
racheal intubation. The swine were then maintained in a deep 
plane of anesthesia with inhaled isoflorane 0.5%-4%. Posi- 
15 tive-pressure mechanical ventilation at a rate of 10-15 
breaths/min and a tidal volume of 500 ml was employed. 

Physiologic transducers were placed as follows. 7.5 French 
sheath introducers (Arrow International, Reading, Pa.) were 
placed in the bilateral femoral arteries. A micromanometer- 
20 tipped catheter (SPC 350, Millar Instruments, Houston, Tex.) 
was fed retrograde to the thoracic aorta from the femoral 
artery for central ABP. The catheter was specifically posi- 
tioned to achieve a diastolic decay that appeared as exponen- 
tial as possible. The second introducer was attached to stiff 
25 fluid-filled tubing (Arrow International) and an external pres - 
sure transducer (TSD104A, Biopac Systems, Santa Barbara, 
Calif.) for femoral ABP. The chest was opened with a midline 
sternotomy. An ultrasonic flow probe was placed around the 
aortic root for gold standard CO (T206 with A-series probes, 
30 Transonic Systems, Ithaca, N.Y.). Finally, a 23- or 25-gauge 
angiocatheter was placed as distal as possible to the brachial 
artery and attached to an external pressure transducer via 
short, rigid tubing for “radial” ABP. Each transducer output 
was interfaced to a microcomputer via an A/D conversion 
35 system (MP150WSW, Biopac Systems). The data were 
recorded at a sampling rate of 250 Hz and 1 6-bit resolution. 

In each animal, a subset of the following interventions was 
performed over the course of 75 to 150 minutes to vary CO 
and other hemodynamic parameters: infusions of volume, 
40 phenylephrine, dobutamine, isuprel, esmolol, nitroglycerine, 
and progressive hemorrhage. Several infusion rates were 
implemented followed by brief recovery periods. 

The technique was applied off-line to six-minute intervals 
(overlapping by three minutes) of the digitized “radial” and 
45 femoral ABP signals to estimate x and proportional CO 
trends. The corresponding gold standard CO trends were 
established by averaging the aortic flow probe measurements 
over the identical time intervals. Gold standard T trends were 
similarly sought by applying the technique of Bourgeois et al 
50 [4] to the central ABP signals. As a metric for comparison 
between an estimated trend (X(i)) and the corresponding gold 
standard trend (X 0 (i)), the root-mean-square-normalized-er- 
ror (RMSNE) in percent was computed as follows: 

55 



60 

where N represents the number of analyzed six-minute inter- 
vals and the argument i denotes the \ th analyzed six -minute 
interval. In order to use this metric to compare an estimated 
65 proportional CO trend with the corresponding absolute gold 
standard CO trend, the former trend was first scaled to have 
the same mean value as the latter trend for each animal. The 
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correlation coefficient (p) between the estimated and gold 
standard trends was also calculated as a scale-invariant metric 
for comparison. 

FIG. 6 illustrates example excerpts of the digitized central, 
“radial”, and femoral ABP signals. The central ABP excerpt is 
from the beginning of the recording period in which the aortic 
catheter was positioned such that the diastolic decay appeared 
as exponential as possible. Unfortunately, the diastolic decay 
of the central ABP signal did not consistently appear as a 
single exponential throughout the recording period in which 
various interventions were employed. The “radial” and femo- 
ral ABP excerpts are from the same time but later in the 
recording period. Although the diastolic decay of the femoral 
ABP excerpt appears smooth, it cannot be adequately repre- 
sented by a single exponential function. 

Table 1 summarizes the results for each animal, and FIGS. 
7 and 8 are examples of the corresponding trends for animal 
1 (worst result) and animal 5 (best result). (The femoral ABP 
signal of animal 4 was not analyzed due to excessive damp- 
ing.) The table and figures illustrate the wide physiologic 
range imposed by the interventions and strong agreement 
between the estimated and gold standard CO trends (in terms 
of CO RMSNE and visually). This strong agreement was 
confirmed by high overall r values (mean±SD) between the 
gold standard CO trends and the CO trends estimated from the 
“radial” (0.84±0.07) and femoral (0.86±0.05) ABP signals. 
Additionally, Table 2 shows that the CO errors (difference 
between calibrated, estimated CO trends and corresponding 
gold standard CO trends) were largely uncorrelated with CO, 
mean ABP (MAP), and heart rate (HR). 
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The low overall x D RMSNE in the table is buttressed by a high 
overall p value (mean±SD) between the two x D trend esti- 
mates (0.85 ±0.08). Thus, there is solid agreement between 
the two x D trend estimates despite substantial differences in 
5 short-time scale morphology between the “radial” and femo- 
ral ABP signals (FIG. 6), and the overall CO measurement 
error is only 14.6%. Preliminary studies show that the six- 
minute intervals used here may be reduced to, for example, 
three-minute intervals without materially affecting the 
results. 

In certain other embodiments of the invention it may be 
desirable to adjust the proportionality constant AC to vary 
depending on the arterial blood pressure or heart rate, if it is 
15 found that there is a systematic difference between depen- 
dence of the estimated cardiac output and the absolute cardiac 
output on heart rate or arterial blood pressure. 

The foregoing description is to be understood as being 
representative only and is not intended to be limiting. Alter- 
20 native systems and techniques for implementing the methods 
of the invention will be apparent to one of skill in the art and 
are intended to be included within the accompanying claims. 
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TABLE 1 


Summary of results. CO is cardiac output; MAP, mean arterial pressure; 
HR, heart rate; RMSNE, root-mean-square-normalized-error; 
and t n long time constant of arterial tree. 


CO RANGE 

MAP RANGE 

HR RANGE _ 

CO RMSNE r%l Tn 

RMSNE 

ANIMAL 

[L/MIN] 

[MMHG] 

[BPM] 

femoral “radial” 

[%] 

1 

1. 6-5.2 

29-100 

96-180 

19.9 

19.1 

8.3 

2 

2. 3-4.2 

54-127 

101-204 

10.2 

16.0 

5.9 

3 

1.9-5. 8 

70-120 

96-186 

8.8 

16.7 

14.7 

4 

1. 3-4.3 

27-106 

103-198 

— 

12.3 

— 

5 

2.4-5. 0 

65-118 

91-198 

10.2 

8.0 

9.7 

6 

2.3-5. 6 

52-108 

109-177 

17.1 

14.7 

12.4 

TOTAL 

1.3-5. 8 

27-127 

91-204 

14.0 

15.0 

10.3 
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Correlation coefficient (p) matrix. ABP is arterial 
blood pressure. See Table 1 caption. 

CO 

MAP 

HR 

“RADIAL” ABP CO ERRORS 0.11 

FEMORAL ABP CO ERRORS 0.09 

0.47 

0.32 

-0.08 55 
0.14 

Since we were not able to obtain gold standard t 

D trends 


from the central ABP signals, we compared the two x D trends 
estimated from the “radial” and femoral ABP signals, which 
should, in theory, be equivalent. Table 1 summarizes the 
comparison results in terms of x D RMSNE in which the gold 
standard trends were established as the average of the two x D 
trend estimates. That is, the x D RMSNE here equals the root- 65 
mean- square of the difference between the two x D trend esti- 
mates divided by their sum and is the same for either signal. 
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We claim: 35 

1. An apparatus for determining cardiac output to within a 
constant scale factor comprising a computer system that 
includes: 

(a) computer readable media having code comprising com- 
puter-executable process steps; and 40 

(b) a processor that executes the process steps to: 

(i) accept an input representing a measurement of an 
arterial blood pressure signal over a plurality of car- 
diac cycles; 

(ii) estimate a function that represents the response of 45 
the arterial blood pressure to a single cardiac contrac- 
tion; 

(iii) fit the function of step (ii) to an exponential-like 
function over a time period that begins a selected 
amount of time following the maximum value of the 
function; 
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(iv) estimate the time constant of the function of step (ii) 
as the time constant of the exponential function of step 
(iii); and 

(v) determine cardiac output to within a constant scale 
factor by dividing average arterial blood pressure by 
the time constant obtained in step (iv). 

2. The apparatus of claim 1, further comprising an analog- 
to -digital converter. 

3 . The apparatus of claim 1, wherein the apparatus includes 
a buffer system. 

4. The apparatus of claim 1, wherein the apparatus includes 
a display device. 

5. The apparatus of claim 1 wherein arterial blood pressure 
is measured invasively or non-invasively at any site in the 
systemic or pulmonary arterial tree. 

6. The apparatus of claim 1 wherein arterial blood pressure 
is further defined as systemic arterial pressure. 

7. The apparatus of claim 1 wherein arterial blood pressure 
is further defined as pulmonary blood pressure. 

8. An apparatus for determining cardiac output to within a 
scale factor comprising a computer system that includes: 

(g) computer readable media having code comprising com- 
puter-executable process steps; and 

(h) a processor that executes the process steps to: 

i. accept an input representing a measurement of an 
arterial blood pressure signal over a plurality of car- 
diac cycles; 

ii. capture long time scale information by estimating an 
impulse response which when convolved with cardiac 
contractions fits the arterial blood pressure signal; 

iii. obtain a time constant by fitting an exponential-like 
function to the estimated impulse response over a time 
period that begins a selected amount of time following 
its maximum value; 

iv. determine cardiac output to within a scale factor by 
dividing a measure of the average arterial blood pressure 
by the time constant. 

9. The apparatus of claim 8 further comprising an analog- 
to-digital converter. 

10. The apparatus of claim 8 wherein the apparatus 
includes a buffer system. 

11. The apparatus of claim 8 wherein the apparatus 
includes a display device. 

12. The apparatus of claim 8 wherein arterial blood pres- 
sure is measured invasively or non-invasively at any site in the 
systemic or pulmonary arterial tree. 

13. The apparatus of claim 8 wherein arterial blood pres- 
sure is further defined as systemic arterial pressure. 

14. The apparatus of claim 8 wherein arterial blood pres- 
sure is further defined as pulmonary blood pressure. 





